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c/3 ■ We argue that in the quarter-filled ladder compound NaV205 the quasi- 

one-dimensional spin system is strongly coupled to a low-energy antiferroelec- 

tric mode of the excitonic type. This mode results from the interplay between 

the electron hopping along the rungs of the vanadium two-leg ladders and 

the Coulomb repulsion between electrons. The charge ordering observed in 
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Ti 
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O ' sodium vanadate at T^ = 34K corresponds to the softening of the antiferro- 

■ electric mode. We consider the spin-isospin model, which describes the spin 
and low-energy charge degrees of freedom in sodium vanadate. Within this 

^ ! model we explain the observed anomalous temperature-dependence of the 

"^ I dielectric susceptibility at T^ and the midinfrared absorption continuum. We 

—I. I identify the broad structure in the low-energy optical absorption spectrum 

0\ • of NaV205 with the continuum formed by two spinous and one low-energy 

^ ■ charge excitation. 
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1 Introduction 

Sodium vanadate (NaV205), first studied in the seventies [|l|, has become 
an object of intensive experimental and theoretical investigations after the 
discovery of the phase transition at Tc = 34K p[. Several physical concepts, 
such as bipolarons , spin-Peierls transition ||2|, |^ , and charge ordering [^, H, 
0, have been proposed to explain the properties of this material. Yet, despite 
all the efforts, the low-temperature crystal structure of sodium vanadate 
remains controversial and the nature of the low-energy excitations in this 
material is still poorly understood. 

In a'-NaV205 the vanadium ions form two-leg ladders organized in layers 
(see Fig. |l|). Until recently it was assumed that the vanadium ladders with 
one electron per rung are equivalent to spin-i chains, so that sodium vana- 
date is a quasi-one-dimensional spin system. Indeed, this material is a good 
insulator and its magnetic susceptibility has temperature behavior similar to 
that of the Heisenberg spin-| chain with the exchange constant J ~ 560K 0. 
At Tc = 34K sodium vanadate undergoes a phase transition into a state with 
a spin gap. Below Tc the lattice period along the ladders (6-direction) and 
perpendicular to ladders within the layers (a-direction) doubles, while the 
period in the direction perpendicular to the layers (c-direction) increases by 
a factor of four 0]. The transition was initially interpreted as a spin-Peierls 
transition, driven by the instability of spin-| chains against the formation of 
local dimers. 

However, more recent experimental studies of NaV205 revealed some in- 
consistencies in that interpretation. First, sodium vanadate does not show 
the strong suppression of Tc by magnetic field, characteristic for spin-Peierls 
systems P, ITDI. Second, the entropy released at the transition is considerably 



higher than the expected release of the spin entropy |||, H, |T0]- Third, the 
dielectric susceptibility, measured by microwave absorption, has an anomaly 
at the transition temperature, which was not observed in the inorganic spin- 
Peierls material CuGeOs [P|. 



Recent X-ray |T^, |1^, |I6[, as well as NMR measurements 0, show that 
above Tc all vanadium ions are equivalent (V^'^"^), whereas below Tc, a charge 
disproportionation occurs. There is, however, still a controversy about the 
number of distinct V sites and their arrangement in the low-temperature 
phase [§|13, 0, ig. 



It was then suggested that the driving force of this transition is the charge 
ordering of electrons in the quarter- filled vanadium ladders and that the spin- 
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gap opens due to alternation of the spin-exchange constants, resulting from 
the charge ordering p, |^. Since sodium vanadate is an insulator both above 
and below Tc, the charge ordering in this material does not originate from the 
Fermi surface instability, as in charge-density-wave materials. Therefore in 
Ref . [0] we proposed a spin-isospin model, in which the rungs of the vanadium 
ladders are assumed to be predominantly occupied by one electron. The 
isospins describe the low-energy charge degrees of freedom, which in sodium 
vanadate are electric dipoles located on the rungs of the vanadium ladders. 
The isospin Hamiltonian (Ising model in a transverse field) is identical to the 
Hamiltonian used to describe (anti)ferroelectric transitions [^. Below Tc the 
dipoles become ordered. Their directions alternate along the ladders [|], 0, 
which corresponding to a zigzag ordering of charges. The antiferroelectric 
nature of charge displacements in the vanadium-oxygen layers was confirmed 
by the observation of anomalies in the dielectric susceptibility close to Tc 



ri| , |T2|, |T^ . In Ref. we showed that the antiferroelectric ordering in the 
system of vanadium ladders opens a spin gap. 

We recall that, depending on whether the variable associated with the 
dipole moment of a unit cell is quantum or classical, there are two limiting 
types of (anti)ferroelectric transitions: a displacive transition and an order- 
disorder transition |2^. If the tunneling between the two states with different 
orientations of the dipole moment can be neglected, the transition is of the 
order-disorder type, described by a classical Ising model. If, on the other 
hand, the hopping between the two states is non-negligible, the transition is 
described by the quantum Ising model in a transverse field. The quantum 
dynamics of the electric dipoles gives rise to a branch of excitations that 
correspond to dipole flips propagating from one unit cell to another. At the 
transition temperature these excitations become gapless. The softening of 
the isospin excitation occurs at the wave vector of the superlattice structure 
appearing in the ordered phase. 

In sodium vanadate the amplitude of the electron hopping along the rungs 
of the vanadium ladders, which plays the role of the tunneling amplitude 
between two states of the electric dipole, is of the same order as the Coulomb 
interaction between electrons on neighboring rungs (leV) that is responsible 
for the charge ordering. Thus the dynamics of the antiferroelectric mode in 
sodium vanadate is essentially quantum and, at least close to the transition 
temperature, there should be low-energy charge excitations in this material. 

For the antiferroelectric ordering (see Sec. ^) the soft mode has a nonzero 
wave vector and, therefore, cannot be directly excited in optical absorption 



and Raman scattering above the transition temperature. Nevertheless, we 
show in this paper that the optical data do provide an evidence for the 
presence of the low-energy charge excitations in sodium vanadate. 

The optical spectrum of sodium vanadate contains the broad absorption 
band covering almost the whole midinfrared region of frequencies |^ , |22| . 
A broad peak was observed also in Raman experiments [^2|, 23, 23]. It was 



suggested in Ref. [^], that the low-energy absorption is due to the pho- 



toexcitation of the two-spinon continuum. However, to convert a photon 
into a pair of spinous, the presence of a permanent electric dipole moment 
is required (both above and below Tc, since the continuum is observed at 
all temperatures). This would be possible if sodium vanadate would have 
the ferroelectric chain-like structure that was originally proposed for this 
material ^, but which later was found to be inconsistent with X-ray data 
I|,|T3|,|T§. 



In this paper we interpret the broad peak in the low-frequency optical 
absorption spectrum as a three-particle continuum: two spinous plus the 
low-energy antiferroelectric excitation. Our mechanism crucially depends 
both on the presence of the soft mode, associated with the charge ordering, 
and on the coupling of this mode to the spin excitations. 

This coupling naturally arises in the spin-isospin model discussed in Sec. |. 
In Sec. ^we consider the charge ordering in this model, which, as was already 
mentioned, is related to the softening of the isospin excitations discussed in 
Sec. ^. In Sec. ^ we apply our model to explain the anomaly of the dielectric 
susceptibility at Tc, while in Sec. || we show how the coupling between spin 
and isospin excitations gives rise to the low-energy optical absorption. We 
discuss our results in the concluding Sec. |^. 

2 Spin-isospin model 

The experimental observations mentioned in Sec. |I] indicate that the low- 
temperature properties of sodium vanadate cannot be described by only con- 
sidering the spin and lattice dynamics, but that it is necessary to include also 
the charge degrees of freedom. On the other hand, sodium vanadate is known 
to be a good insulator, both above and below Tc. We, therefore, argued in 
that the relevant electronic excitations in this material are excitons. In 
this section we briefly discuss the "spin-isospin" model, introduced in Ref. 0, 
which describes the low-energy excitations in sodium vanadate. 




Figure 1: The crystal structure of the V-0 plane in NaV205: oxygens are located at the 
corners of the plaquettes, while V ions are located at their centers; the shaded plaquettes 
are vacant. Also shown are the vanadium ladders (dashed lines) and the relevant dxy~ 
orbitals of the V ions. 



In that model we restrict ourselves to the states with only one electron 
per rung of the vanadium ladder. This was first suggested by Smolinski et al. 
on the basis of LDA calculations [0 , which show that the electron hopping 
amplitude along rungs, t^, is significantly higher than the amplitudes of 
hopping between the sites of different rungs (of which the largest is the 
hopping amplitude along ladders, t\\). If one first neglects the interrung 
hopping and the interrung electron-electron interactions, then on each rung 
of the quarter-filled vanadium ladder there is one electron occupying the 
symmetric bonding state with energy —t±. Due to the large value of the 
on-site Coulomb repulsion (which we assume to be infinite), an electron can 
hop only to an empty site of a neighboring rung. Then the hopping to a 
neighboring rung results in an energy increase 2t^. Thus, the quarter- filled 
two-leg ladder becomes equivalent to the half-filled Hubbard chain with an 
effective "on-rung" repulsion Ur = 2t±. 

This explains, in principle, why vanadium ladders are insulating. How- 
ever, for the values of the hopping amplitudes, obtained in [|l^, t± = 0.38eV 
and t|| = 0.17eV, the width of the one-dimensional band, 4t\\, is comparable 
with Ur- Thus the charge fluctuations on the ladder rungs are not small and 



the treatment of the vanadium ladder as a spin chain, as it was done in 
is not well-justified. 

The situation improves if we take into account the Coulomb interactions 



between electrons on different rungs. The latter lead to an increase of the 
value of the "on-rung" repulsion. For instance, if we include the repulsion 
between the nearest-neighbor sites in ladders, V, then Ur ~ 2t± + ^, for small 
V. In general, the interrung Coulomb interactions make Ur not a well-defined 
quantity, as it becomes dependent on the positions of many electrons. 

Thus, on the one hand, the interrung Coulomb interactions help to justify 
the assumption of one electron per rung. On the other hand, they strongly 
mix the symmetric and antisymmetric states on each rung, as the typical 
value of such interactions is of the same order as the energy separation be- 
tween these two states, 2t±_ ^ 0.75eV. Thus, the electron position on a rung 
becomes an important additional degree of freedom, which makes the vana- 
dium ladder different from a spin chain. 

The two states, «„ and dn, corresponding to the two possible positions 
of a single electron on the rung n, can be described as the up and down 
eigenstates of an isospin-^ operator, T^. With two spin projections there 
are in total four different states of an electron on a rung (see Fig. ^. The 
spin Sn and the isospin Tn are defined on the lattice, the sites of which 
correspond to the rungs of vanadium ladders (see Fig. ^. In this paper we 
ignore the three-dimensional structure of sodium vanadate and only consider 
one oxygen-vanadium layer. Then the lattice of our model is triangular. 
Furthermore, for our considerations it will not be important that in each layer 
the orientations of pyramids in the nearest-neighbor ladders are opposite 
to each other (which is why the unit cell above T^ contains two ladders). 
Therefore, we chose f i = b and f2 = |(a + b) as the basis of the unit cell in 
the high-temperature phase. 

The charge part of the Hamiltonian of the model, which only includes the 
isospin operators, has the form: 

Ht = -2t^ E ^n + ^ E Vnm ^ T^, (l) 

n ^ nm 

where t±_ is the amplitude of electron hopping along the ladder rungs and the 
second term describes the Coulomb repulsion between electrons on different 
rungs. The amplitudes Kim are defined by 

Km = K(nn, Urn) + K(dn, d^) - K(nn, d^) - K(dn, Urn), (2) 

where, e.g., V{un, dm) is the energy of the Coulomb interaction between an 
electron on rung n in the up state and an electron on rung m in the down 
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Figure 2: The four different states of a single electron on a rung represented as the 
eigenstates of the spin S — ^ and the isospin T = ^ operators. 
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Figure 3: From vanadium ladders (left-hand side of the picture) to the effective lattice 
(right-hand side of the picture) . The sites of the latter should be identified with the centers 
of the rungs of the V ladders. 



state. The Coulomb interactions favor a certain charge ordering, which in 
our model corresponds to an ordering of the ^-projections of the isospins, 
while the electron hopping counteracts the ordering. Since the hopping term 
does not commute with the Coulomb term, the isospins are quantum (rather 
than classical) degrees of freedom, like the orbital degrees of freedom in Mott 



insulators with orbital degeneracy p5 



We now turn to the spin contributions to the Hamiltonian. The spin 
exchange results from virtual hopping of electrons between neighboring rungs, 
and the exchange interaction is the strongest in the ladder direction. For the 
half-filled Hubbard chain the exchange interactions between spins on two 
neighboring sites has the form jj-Si2, where S12 = 2 (Si ■ S2) + | is the 
operator that exchanges Si with S2. In the two-leg ladder it is possible 
to have in the intermediate state a doubly occupied rung without doubly 
occupied sites. Thus even for infinite U the exchange coupling is nonzero and 



U has to be replaced by the "on-rung" repulsion Ur- In the infinite U case 
the exchange of spins necessarily involves the exchange of the 2;-projections 
of the isospins on two neighboring rungs {uid2 <-*■ diU2), so the spin-exchange 
interaction has the form: ^Su (T^T2 + T^T2~]. Thus, instead of a pure 
spin-exchange Hamiltonian we obtain a spin-isospin interaction: 

^S = A^(S„-S„+fJ (T„^T„^+f,+T„^T„\fJ . (3) 

n 

Here, the rungs n and n + fi are the two nearest-neighbor rungs of the ladder 

QJ-2 

(see Fig. ^) and the interaction amplitude A = jj^, Ur being the effective 
"on-rung" Coulomb repulsion discussed above. 

We note that Eq. (^ gives only an approximate description of the spin- 
isospin interaction in sodium vanadate, because in the presence of inter- 
rung Coulomb interactions, both initial and intermediate states in the spin- 
exchange process are complicated many-electron states, so that the increase 
of the energy in the virtual state cannot be described by a single constant Ur- 
However Eq.(|D contains many basic ingredients, necessary for a qualitative 
description of the properties of sodium vanadate. 

In particular, treating the spin-isospin interaction in the mean-field ap- 
proximation we obtain the effective spin-exchange coupling constant: 

J = A(T„-T„\f^+T„^T„\f^), (4) 

where the brackets denote the thermal and quantum average of the isospin op- 
erators. The exchange coupling is, therefore, temperature-dependent, which 
may account for the deviation of the magnetic susceptibility of sodium vana- 
date from the Bonner-Fisher curve |jlO|, |26|. Using the Hamiltonian Eq.(|l]), 



one can show that, above Tc, the exchange coupling increases as the temper- 



ature decreases, which is in agreement with the conclusion of Ref. I^g, where 
such an increase was introduced phenomenologically in order to reconcile the 
magnetic susceptibihty data for NaV205 with the temperature dependence 
of the magnetic susceptibility of the Heisenberg spin-| chain. 

However, there is one very important effect of the electronic charge dis- 
tribution on the spin exchange that is not described by Eq.(^. Namely, the 
Hamiltonian Eq.(^ does not explain the spin-gap opening due to the zigzag 
charge ordering. In order to find the term in the spin-isospin interaction re- 
sponsible for the spin-gap opening, we note that the exchange between spins 
on rungs n and n + f i depends on the position of electrons in the nearest 



rungs of the two neighboring ladders, n + f 2 and n + f 1 — f 2 (see Fig. (j^)). 
Roughly speaking, electrons affect the hopping amplitudes between the rungs 
of nearest-neighbor ladders, which in turn, affect the spin-exchange. The cor- 
responding Hamiltonian describing such an influence has the form: 

HsT = ^ H (^n+f2 ~ ^n+fi -f2 ) (Sn ■ ^n+fi ) • (5) 

n 

As was argued in Ref. 0, this interaction is responsible for the opening of 
the spin gap in the charge ordered phase. Indeed, for the zigzag structure 
with a doubling of the period in the a and h directions (see Fig. ^a), the 
exchange alternates along ladders, which opens a spin gap. 

To end this section, we note that the virtual hopping of electrons on 
neighboring rungs, apart from affecting the spin-exchange, also results in 
corrections to the pure isospin Hamiltonian Eq.(|l]), which we neglect here, as 
they are relatively small. Furthermore, for large, but finite on-site Hubbard 
U , one obtains also the spin-isospin interaction term of the form: 

^l? = ^'E(Sn-S„+fjT„^T„V^, (6) 

n 

where A' oc tn/f/. 

3 Charge ordering 

In this section we ignore spins and discuss the phase diagram for the Hamilto- 
nian (|l]), which describes the charge dynamics. This is justified as, according 
to our assumption, the driving force of the transition in sodium vanadate is 
the electronic charge ordering. 

Since the Hamiltonian (|I|) is real, the average value of T^ is zero at all 
temperatures (unless the time-reversal symmetry is spontaneously broken in 
the ordered state). On the other hand, M^ = (T^) 7^ at all temperatures, 
due to the presence of the "external" transverse field 2t±. The ordered state 
of the isospin system is characterized by a nonzero value of M^ = (T^). The 
ordering of isospins in the z-direction corresponds to the modulation of the 
electronic density on the rungs of the vanadium ladders. The type of charge 
ordering depends on the details of the Coulomb interactions Kim- 

If we use as estimates the unscreened Coulomb interactions between elec- 
trons, then, according to Eq.(^, the Ising interaction between the isospins 
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Figure 4: Two charge ordered states: zigzag (or antiferroelectric) ordering with the 
doubling of the unit cell in the directions a and b (a) and the chain-like (ferroelectric) 
structure (b). The left part of each picture shows the distribution of the electronic density, 
while the right part shows the corresponding ordering of isospins. 



on the nearest-neighbor rungs in the ladder, Ja = ;^Ki,n+fi = 0.62eV. This 
"antiferromagnetic" interaction favors a zigzag structure, hke e.g. the one 
shown in Fig. ^. On the other hand, the interaction between the isospins 
on the nearest rungs of two neighboring ladders, Jp = ;|Ki,n+f2 = 0.54eV, is 
"ferromagnetic" (for simplicity, we neglect the shifts of vanadium ions out of 
the basal plane). The "ferromagnetic" interaction favors the structure shown 
in Fig. ^, in which electrons occupy with higher probability one chain of each 
ladder. In general, the charge ordering corresponds to the ordering of the 
isospins in the 2;-direction: 

(T^) = e*^-^"M", (7) 

where Xn is the position of rung n, and Q = (-p, -p) for the zigzag structure, 
while for the chain structure Q = 0. 

If one would take into account only the Ising interactions between the 
neighboring rungs, then for .J a > \Jf\ the energy of the zigzag structure 
is lower energy than that of the chain structure. However, since the differ- 
ence between J a and \Jf\ is relatively small, the effects of screening of the 
Coulomb interactions and the long-range interactions between the isospins 
are important for the determination of the nature of the charge-ordered state. 

Though the precise calculation of the Ising interactions between the iso- 
spins is very difficult, one can discriminate between different charge-ordered 
structures on the basis of the available experimental data. In particular, the 
chain-like ordering (see Fig. ^d) does not result in the increase of the lattice 
periodicity in the a and b directions observed in NaV205 below T^. Further- 
more, the chain structure corresponds to a ferroelectric state. As we shall 
show below, the optical response of such a system would be quite different 
from that of sodium vanadate. On the other hand, the zigzag charge order- 
ing results in a doubling of the lattice period in the 6-direction. Moreover, if 
the zigzags in the next-nearest-neighbor ladders have opposite phases (as in 
Fig. Ha), then the lattice period in the a-direction also doubles. 

The generic phase diagram of the Ising model in a transverse field is shown 
in Fig. |[ For large values of t± the system is disordered at all temperatures. 
For t± < t*_i, the system becomes ordered below some critical temperature 
Tc. At t± = t*j_ the system has a quantum critical point, at which Tc = 0. 
Since in sodium vanadate the typical values of both the Ising couplings and 
the hopping amplitude t± are of the order of 0.5eV, whereas T^ is only 34K, 
this material has to be very close to the quantum critical point. 
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Figure 5: Generic phase diagram of the model described by Eq.(P in the T-t±^ plane. 

The proximity to the quantum critical point may be responsible for the 
strong dependence of Tc on the hydrostatic pressure observed in sodium 
vanadate p7| , |25[] . We note that the fact that the transition temperature 
is suppressed by pressure can be naturally explained in our model. Both 
the hopping amplitude t± and the amplitudes of the dipolar interactions be- 
tween the rungs, Kim, increase under pressure. However, t±, which involves 
the overlap of electronic wave functions, depends on the distances between 
ions more strongly then Vnm- Since the hopping amplitude tends to disorder 
the system, the transition temperature decreases under pressure. 

The isospin Hamiltonian Eq.(|I|) is the Hamiltonian of the Ising model with 
long-range interactions in a transverse field. Recently, a considerable progress 
was made in the calculation of the linear response functions of such systems, 
which close to a quantum critical point, become to a large extent p9[. While 



those results could be useful to describe the response function at wave vectors 
q close to Q {e.g., the temperature dependence of the shape of neutron 
diffraction peaks), in this paper we are mostly interested in the dielectric 
and magnetic susceptibility of sodium vanadate at q = 0, as the latter is 
measured in optical experiments. We use the random phase approximation 
(RPA) to obtain the necessary correlation functions of the isospin operators. 
Then, for our purposes, it is sufficient to describe the ordered state in the 
mean field approximation. 



11 



In this approximation M^ and M^ are found from 
2t^M^ + M^ Em KmM^ = 0, 



(8) 



M„ = J{M-y + {M^Y = \ tanh ^. 



The first equation fixes the orientation of the vector (Tn), which has to 
be parallel to the vector of the "field" applied to the isospin at the site n, 
Hn = (2tx, 0, — Em ^nm^m)- The secoud equation gives the length of (Tn) 
(in that equation H^ = |Hn|). 
Above Tc we have 

r M^ = ifanh^, 

(9) 

M^ = 0, 

while for the ordered states shown in Fig. ^, the order parameter is given by 
Eq.(^ and M^ is constant below T^. 

M"" = w = — ^. do) 

|\/(Q)| ^ ' 



The absolute value of the magnetization, M = J{M^Y + (iVf^)^, satisfies 



2 2T ^ ' 

where ^(Q) is the value of the Fourier transform of V^m at q = Q: 

\/(Q) = 5:e-*Q(^"-^-V„m. (12) 

n 

Finally, at the critical temperature 

2tx + M^X^c) V^(Q)=0, (13) 

where M^{T) is given by Eq.(|^). Clearly, V^(Q) has to be negative and ^(q) 
has to have a minimum at q = Q (the location of the minimum determines 
the type of the charge-ordered structure). Furthermore, Eq.(0) only has a 
solution for |V^(Q)| > At^. In other words, in the mean field approximation, 
the critical point is reached at t\_ = ■'— x^. 
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4 Isospin excitations 

The isospin excitations described by the Haniiltonian (Q) are the rung exci- 
tons, i. e. the electronic excitations from the symmetric state with energy —t± 
to the antisymmetric state with energy +t±, which can propagate from rung 
to rung due to the Coulomb interactions between the electrons. This can be 
most easily understood by performing a rotation in the isospin space around 
the y-axis, under which T^ —>■ T^ and T^ — >■ —T^ (see Appendix^. Then the 
excitation from the symmetric to the antisymmetric electron state on a rung 
corresponds to a spin flip. Due to the Ising interaction between the isospins 
on different rungs, which in the new basis has the form, | I]nm Kim ^n^mi 
the spin flip can hop from rung to rung. 

The dispersion of these excitations can be found, e.g., by considering the 
retarded Green function 

POD 

((T"|T"))^q = -i / c/te*("+*^)*y e-*'^-"" {\T'{t),T\0)]) (14) 

Jo V 

(here we use the original basis in the isospace). This Green function is 
calculated in the random phase approximation (RPA) in Appendix ^ (see 
Eqs. ( p3D and (1^)). From the poles of the RPA Green function we obtain 
the energy of the isospin excitation with wave vector q, valid both above and 
below Tf.: 



^q = \/{2t±y + 2t^M-V{q) + {\V{Q)\M^y, (15) 

where the averages M^ and M^ are given by Eq.(|^) above T^ and Eqs.(p!OD 
and (|TTp below Tc. 



Since the energy of the electron excitation on a rung, 2t± ~ 0.7eV |T6 
is rather large, the isospin excitations are irrelevant for the low-temperature 
properties of sodium vanadate, unless this large energy is compensated by a 
large band width (the latter is determined by the amplitude of the variation 
of V{q), which for sodium vanadate is also expected to be of the order of 
leV). We thus assume, that close to the bottom of the exciton band {i.e., for 
q ~ Q), the exciton energy is comparable to temperature (or smaller). 

The existence of soft isospin excitations implies a strong temperature de- 
pendence of the exciton gap, A = Eq. The typical temperature dependence 
of the isospin gap is shown in Fig. ^. Above Tc the "transverse magnetiza- 
tion" M^ grows when temperature decreases (see Eq.(P)), which results in 
a decrease of the gap (we recall that K(Q) in Eq.(|l5|) is negative). At the 
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Figure 6: Temperature dependence of the isospin gap A. The value of |y(Q)|, used 
for this calculation, is close to its quantum critical value, 2t±, so that the isospin gap 
A<2ij.. 



transition temperature, Tc, the exciton gap disappears in accordance with 
Eq.(|l^). Below T^ the order parameter, M^, increases as temperature de- 



creases, while M^ stays constant (see Eqs. ( [T0| ) and (|TT])). As a result the 
exciton gap again reappears. 



5 Anomaly of dielectric susceptibility at Tc 

The isospin excitations discussed in the previous section can be excited by 
applying an electric field in the direction of the rungs (a-direction). Within 
the isospin model the interaction with such a field, Ea{t), is described by the 
Hamiltonian 

H,^t = elEait)J2T^, (16) 



where / ~ | is the length of the rung. Then the contribution to the dielectric 
susceptibility of the electrons occupying the dxy orbitals on the vanadium 
sites is: 

Xa{co,cD = -C{{T^\T^))^,^, (17) 

where the constant C = ^^^, a, b, and c being the lattice constants of 
the high-temperature unit cell, and the retarded Green function ((T^|T^))^_q 
above and below Tc is given by Eqs. ( ^3|) and (42), respectively. Using these 
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equations we obtain 

Xa(u;, q) = -C "^ ^ (18) 

for the charge susceptibihty above Tc (here, M^ is given by Eq.(||)) and 

below Tc- 

The frequency dependence of the dielectric susceptibility was measured by 
optical absorption ||2T|, ^ |5D[ . Its imaginary part shows a rather broad peak 
at energy ~ leV, which was assigned in Ref. |21|] to bonding-antibonding 



excitations of electrons on rungs. According to Eg . ([181) the peak should be 
located at the energy of the isospin excitations with q = 0. As for the 
zigzag structure, shown in Fig. 4a, this excitation lies far from the bottom 
of the exciton band, located at Q = (-p, -p), its energy, -Eq=o, should indeed 
be rather large (~ 2t^). On the other hand, for the chain structure (see 
Fig. 4b) the bottom of the exciton band would be at q = and the soft 
isospin excitations would result in a strong peak in the optical absorption at 
low energy, -Eq=o ~ 0, which is not observed in these experiments. 

The type of charge ordering also determines the type of the anomaly in 
the temperature dependence of the static dielectric susceptibility close to Tc- 
Such an anomaly appears due to the redistribution of charge at the phase 
transition. The order parameter M^, as well as M^ and the isospin en- 
ergy Ec^ depend on T, which makes the dielectric susceptibility temperature- 
dependent (see Eqs. ([l8|) and (pjf )). 

For the antiferroelectric ordering of charges, as in the case of the zigzag 
ordering (see Fig. pi), the susceptibility to the uniform field is finite at Tc 
(see Eq.(p!8|) with -Eq=o finite), while its temperature derivative has a dis- 
continuity at the transition temperature, as is shown in Fig. 0. (This is 
analogous to the anomaly in the magnetic susceptibility of the Ising antifer- 
romagnet). On the other hand, for the ferroelectric type of ordering, such as 
the chain-like structure shown on Fig. ^, the susceptibility would diverge at 
Tc- The anomaly in the temperature dependence of the dielectric susceptibil- 
ity of NaV205 has been measured at microwave and far-infrared frequencies 
ri| , p!2| , |T^. The observed anomaly is similar to that shown in Fig. |^ and 
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Figure 7: The typical temperature dependence of the static dielectric susceptibility 
to electric fields applied in the a-direction calculated in the RPA approximation. The 
parameters of the isospin Hamiltonian were chosen such that the value of Tc is much 
smaller than t±. 



consistent with a charge ordering of the antiferroelectric type |Tl| , although 
above Tc the dielectric susceptibility of sodium vanadate continues to grow 
with temperature, whereas the RPA predicts a decrease. 

6 Optical absorption at low energies 

A remarkable feature of the optical absorption spectrum of sodium vanadate 
at low-energies is that the low-energy shoulder of the peak at 8000cm~^ 
extends all the way down to ~ 3000cm^^ (see Fig. ^). Furthermore, at 
lower frequencies there is another broad absorption band stretching from 
~ lOOcm"^ up to ~ 1500 — 2000cm~^ with a maximum at ~ 300cm~^ ||21|, |22 . 
A very broad peak with a maximum at 600cm~^ is also observed in Raman 
scattering ||2^, ^, ^]. The large width of these peaks, as well as a number 
of the Fano resonances observed in the far-infrared absorption spectrum ||2^, 
2^ , indicate the presence of a broad continuum of low-energy excitations in 



sodium vanadate, which covers, practically, the whole midinfrared range of 
frequencies. 

The low-energy continuum, observed both above and below T^, can, in 
principle, be due to the photoexcitation of several isospin excitations lying 
close to the bottom of the exciton band. To describe such processes one has to 
go beyond the RPA. Then the one-exciton states are mixed with the states 
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Figure 8: The frequency dependence of the real part of the optical conductivity of 
NaV205 at 300K in the midinfrared range (courtesy of Prof. D. van der Marel [2l]). 



containing three, five, and more excitons. However, since the number of 
excitons created in the photoabsorption is odd, the momentum conservation 
does not allow all of them simultaneously to have low energy. Thus, the pure 
exciton continuum should begin at rather high energies and cannot account 
for the mid-infrared continuum in sodium vanadate. 

To substantiate the point that the isospin excitations alone cannot give 
rise to the observed low-energy optical absorption, we consider an isolated 
vanadium ladder with the isospin Hamiltonian: 

n n 

For open boundary conditions, this Hamiltonian can be diagonalized exactly. 
To this end we first perform a rotation around the y-axis in the isospin space 
(c/. Appendix). The resulting Hamiltonian is diagonalized by transforming 
the isospins to fermions (using the Jordan- Wigner transformation) and then 
applying the Bogoliubov transformation [^. The exact dispersion of the 
fermionic excitations reads: 

E„ , ___^, , , , 

2 V V2 

with q a quantum number that follows from a transcendental equation |^ 
and which lies in the interval < g < vr. Thus, the isospin gap is given by 

V 




A = E^ 



2tx- 
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Figure 9: The oscillator strength for an isolated two- leg ladder of A^ = 70 rungs with 
very small value of the isospin gap at g = tt (see explanations in the text). The narrow 
peak is due to the single exciton absorption. The absorption above the main peak is 
due to the photoexcitation of three excitons. The very small oscillator strength for the 
frequencies below the frequency of the main peak is due to the one-exciton absorption and 
is a finite-size effect. 



The oscillator strength of a n-exciton state (n odd) may now be expressed 
in terms of the ground state expectation value of strings of Fermi opera- 
tors. The latter may be evaluated using Wick's theorem which results in 
determinants of matrices with components that follow from the Bogoliubov 
transformation coefficients. For the 1- and 3-exciton states, this calculation 
has been performed explicitly in Ref. |^. Using these results, we find for 
a ladder of A^ = 70 rungs and 2t± = 0.7eV the absorption spectrum plot- 
ted in Fig. 1^. As NaV205 seems to be close to the quantum critical point, 
2t± = ^, (see discussion in Sec. 3), we choose ^ = 0.6999eV, which gives 
an isospin gap 10~^eV. The narrow peak is due to the single-exciton absorp- 
tion. The absorption above the main peak is due to the photoexcitation 
of three excitons. This process contributes to the broad high-energy wing 
of the main absorption peak observed in sodium vanadate pl|. The small 



oscillator strength at frequencies below the frequency of the main peak is 
due to the one-exciton absorption and is a finite-size effect (breaking of the 
momentum selection rule). For small photon frequencies the absorption due 
to the 5-exciton states is negligible due to the small phase volume of these 

states. 



In Ref. 1 21] the low-energy optical was associated with the two-spinon 



continuum. Since magnetic excitations cannot be directly induced by photon 
absorption, the authors of Ref. ||21| invoked the so-called "charged magnon" 



mechanism. The corresponding Feynman diagram is shown in Fig.^(a): The 
process involves the creation of a virtual high-energy isospin excitation with 
q = 0, which then decays into two spinous. However, the coupling of the 
isospin excitation with zero wave vector to spinous is only nonzero if the 
system has a permanent electric dipole. In other words, the magnons have 
charge only if sodium vanadate is ferroelectric, which would require the chain- 
like electronic ordering (cf. Fig 4b) both above and below Tc. For a long 
time sodium vanadate was, indeed, believed to have such a chain structure 
[jl|. However, this structure was found to be inconsistent (certainly, in the 
high-temperature phase) with recent experimental data (see discussion in 
Sec. 5). 

Thus the coupling of a photon with zero wave vector to two spinous is 
zero, and we have to find another mechanism for the low-energy optical ab- 
sorption. In the remainder of this section, we shall show that the low-energy 
absorption may result from the photoexcitation of a three-particle contin- 
uum: two spinous and one low-energy isospin excitation. This mechanism is 
illustrated by the diagram in Fig.[Ty(b). This process also takes place via a 
virtual isospin excitation with zero wave vector and high energy Eq, which 
then decays into a low-energy isospin excitation with the wave vector q close 
to the bottom of the band, q ~ Q, and two spin excitations, which carry the 
momentum — q. Alternatively, the thermally excited isospin excitation can 
be annihilated in the process of optical absorption (see Fig.p!0|(c)). 

Our mechanism requires the coupling of two isospin excitations (one of 
high energy, another of low energy) to two spin excitations. In sodium vana- 
date such a coupling results from the fact that the spin exchange of electrons 
on neighboring rungs of a two-leg ladder cannot be separated from the ex- 
change of their isospins. The Hamiltonian describing the spin-isospin cou- 
pling was discussed in Sec. 0. As an exact treatment of the spin-isospin model 
is not possible we use a number of approximations to calculate the optical 
absorption spectrum. First we note that the Hamiltonian of the spin-isospin 
interaction given by Eqs. (^, (^), and (^) can be written in the form: 

UsT — ^ST + ^ST + ^ST — Z^ -'n'->n, 

n 

where 7^ contains all isospin operators: 

r„ = A (T„-T;^Vf. + ^^^n+f J + A'T^T^^,^ + B {t^^,^ - T^^,,^,^) , (20) 
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Figure 10: Diagrams describing: (a) the "charged magnon" mechanism, (b) the pho- 
toexcitation of the two-spinon continuum and a low-energy exciton, and (c) the photoex- 
citation of the two-spinon continuum with the annihilation of a low-energy exciton. The 
wavy and dashed lines correspond to, respectively, a photon and a charge exciton, while 
the black square indicates the two-spinon continuum. 



and 

Sn = bn'^n+fi- 

Although in our model the spin dynamics is inseparable from the dynam- 
ics of the isospins, we may, however, obtain a "pure spin Hamiltonian" by 
treating the isospin operators in Hst within the mean-field approximation: 

Hs = / .'-^n.n+fi (Sn-Sn+fi) , (21) 

n 

where the effective spin-exchange constant is obtained by the thermal average 
of Tn over the state of the isospin system: 
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Figure 11: The diagram describing the correction to the Green function of the isospin 
excitation due to scattering on the spin excitations. This correction which gives rise to 
low-energy optical absorption. 

As was already discussed in Sec. ||, the thus obtained spin-exchange constant 
depends on temperature and becomes alternating along the spin chains below 
Tc. Since according to our assumption the charge ordering is mainly driven 
by the pure isospin Hamiltonian Ht, we neglect the effect of the spin-isospin 
interactions on the state of the isospin system. The remaining part of the 
spin-isospin Hamiltonian (the residual interaction) describes the scattering 
of isospin excitations on spin excitations. This scattering gives rise to the 
low-energy absorption and will be treated as a perturbation. 

In this paper we consider the optical absorption only above T^ {i.e., no 
charge ordering and no spin gap). Furthermore, the isospin excitations are 
treated as bosons. In particular, the operator 

^ y/Nn 

annihilates the isospin excitation with the wave vector q and creates such an 
excitation with the wave vector — q: 



The coefficient in the right-hand side of the last equation was found by com- 



parison of this equation with the Green function Eq. fpSl) . 

Since the operator 7^ (see Eq.(pOD) contains products of two isospin op- 
erators, it either results in the creation and annihilation of two isospin ex- 
citations, or in the scattering of the isospin excitation. In particular, if one 
of the isospin excitations has zero wave vector (as is the case in the optical 
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absorption), then 






where f^q is the couphng constant. 

To lowest order in the residual spin-isospin interaction, the Green function 
of the isospin excitation obtains the correction shown in Fig. |TI|, which leads 
to the low-energy optical absorption described by the diagrams in Figs. |l^ 
(b) and (c). The corresponding contribution to the imaginary part of the 
dielectric susceptibility is: 



^'^"("•'" - f?S 



Xsi^ - Eq, -q) coth -^ - coth "^ 



+ x'si^ + E^,ci)(coth^-coth^^''^'' 



2T 2T J 

(22) 



2T 2T J 

Here, the first term in square brackets is due to photoexcitation of one low- 
energy exciton plus two spinous (see Fig. [IO|b), while the second term de- 
scribes annihilation of a thermally excited low-energy excitation, accompa- 
nied by the creation of two spin excitations (see Fig. p!0|c). In Eq.(p2D we 
have defined 

1 ( 2eltlM^ V 

ahc\ El -u^J ' 
while Xsi^j q) = Xs{^, Qi) is the susceptibility of the Heisenberg spin-| chain: 



oo 



Xsiu;,q,) = I / rfte^-*-^''^"([5„(t),5o(0)])^, (23) 

Jo 

(here gi is the wave vector in the chain direction measured in units of \). 
Since no exact expression for this susceptibility is known that holds for all 
values of g, w, and T, we calculate it by substituting the Heisenberg model 
by the renormalized XY model, which is equivalent to a half-filled chain of 
spinless fermions with the dispersion e{k) = — pj cos k, where p = 1 + 2/7r 
In. Then 

and 
Xsi^^Q) = ^ J dk6 {uj + Sk - Sk+q 
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According to our assumption, the isospin excitations only are soft in a small 
vicinity of gi = vr. Then the wave vector of the spin excitations in the low- 
energy continuum is also close to vr, in which case the analytical expression 
for the imaginary part of the spin susceptibility is: 




u; \^ , uj 



xU^,n) = e{Q'-u')^^ll - y tanh— . (24) 

Here, Q = 2pJ is the bandwidth of the spinless fermions, which corresponds 
to the maximal energy of the two-spinon excitation in the Heisenberg model. 
In Fig. |T^, we show the real part of optical conductivity, ai^u) oc uXai^, 0), 



resulting from Eqs. (0) and (|2^) for the model parameters: A = 200K, 
T = 300K, and J = 560K is shown Fig. ^. The spectrum has two singu- 
larities at a; = fi ± A due to the fact that the maximal energy carried by 
two spin excitations cannot exceed Q. The contribution of the process with 
the annihilation of an isospin excitation (see Fig. p!0|c) has its maximum at 
uj ^ Q and decreases exponentially (as e~^'^) at larger frequencies. There- 
fore, the absorption at fi, T -C cu -C ii^o is entirely due to the process shown 
in Fig. p!0|b. At these frequencies the optical absorption spectrum measures, 
basically, the density of isospin excitations. Since the spectrum of these ex- 
citations stretches from very low energies near the bottom of the band at 
q = Q up to the energies of the order of 2t±, a (relatively small) optical ab- 
sorption due to the photoexcitation of the three-particle continuum should 
occur at all frequencies between A and Eq. This explains the midinfrared 
absorption continuum observed in NaV205 and described in the beginning of 
this section. We note that, in agreement with our theory, this continuum is 
only observed when the electric field is directed along the ladder rungs (along 
a). 

To end this section, we point out, that our mechanism of the optical 
absorption at low photon energies is very similar to the Lorenzana-Sawatzky 
mechanism of the phonon- mediated photoexcitation of two magnons P5| . 
The difference lies in the fact that the phonon dispersion is usually small 
or comparable to that of magnons, so that the shape of the three-particle 
continuum is mainly determined by the spin excitations |3^. On the other 



hand, in our case the shape of the optical spectrum is governed by the isospin 
excitations, which a have much wider band than spinous. 
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Figure 12: The low-energy part of the optical absorption spectrum due to the process 
shown in Figs.^(b) and (c). 

7 Discussion and Conclusions 

In this paper we discussed the spin-isospin model of sodium vanadate. We 
showed that when the rungs of the vanadium two-leg ladders are occupied 
mostly by a single electron, the spectrum of low-energy excitations may in- 
clude, not only spin, but also charge excitations. The latter are the bonding- 
antibonding electronic excitations on the ladder rungs, which become itiner- 
ant due to the Coulomb interactions between electrons on different rungs. 

We described the dynamics of these excitations with isospin-^ operators. 
The isospin Hamiltonian is identical to the Hamiltonian of the Ising model in 
a transverse field. In this model one obtains a transition into an ordered state 
al low temperatures, which corresponds to the charge ordering of electrons on 
rungs of the vanadium ladders. At the transition temperature the spectrum 
of the isospin excitations becomes gapless. 

Furthermore, we showed that the dynamics of the spins and isospins are 
coupled to each other. A pure spin exchange does not exist - the exchange of 
spins necessarily involves the exchange of isospins. As a result the effective 
spin exchange constants depend on temperature. Below the charge order- 
ing temperature the spin exchange constants alternate along ladders, which 
opens a spin gap. 

As we have shown, the spin-isospin Hamiltonian explains the anomaly 
in the dielectric susceptibility at Tc and the optical absorption at low fre- 
quencies, observed in sodium vanadate. The shape of the anomaly in the 
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dielectric constant proves that the charge ordering below T^ in NaV205 is of 
an antiferroelectric (zigzag) type. 

A considerable part of this paper was devoted to the calculation of the 
optical absorption within the spin-isospin model. We argued that neither 
pure isospin excitations, nor pure spin excitations can explain the broad 
bands in the optical absorption spectrum of sodium vanadate found in the 
whole range of mid-infrared frequencies. We showed that this absorption 
can be due to the photoexcitation of a three-particle continuum (one low- 
energy exciton and two spinous). We also note that below Tc, the quasi- 
momentum conservation allows for the optical excitation of a single low- 
energy excitation, as the increase of the lattice period folds these modes into 
the zone center. Indeed, both in optical absorption p2| and Raman spectra 



24| , several lines have been observed whose intensity shows an anomalous 
temperature dependence. 

Thus the spin-isospin model allows us to explain qualitatively the main 
experimental facts on sodium vanadate, i.e., the charge ordering, the spin-gap 
opening, and the low-energy absorption. The main difficulty of this model is 
to explain the low value of the charge ordering temperature. One possibility 
mentioned above is the proximity to the quantum critical point, which, how- 
ever, requires a fine tuning of the model parameters (the transverse field and 
the Ising interaction between the isospins). Another possibility, the frustra- 
tion of the Ising interactions due to the relative shifts between neighboring 
ladders, will be discussed in detail elsewhere. 

In conclusion, in sodium vanadate the electronic charge dynamics seems 
to play an important role, and the coupling between the charge and spin 
degrees of freedom is crucial in understanding the properties of this material. 
Finally, we note that a similar spin-isospin model could be used to describe 
the interplay between the charge and antiferromagnetic ordering in the quasi- 
two-dimensional organic system k-(BEDT-TTF)2X (see, e.g., [^ i8|), in 



which the role of the V rungs is played by dimers of BEDT-TTF molecules. 
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A Spectrum of isospin excitations in RPA 

In this Appendix we find the spectrum of isospin excitations both in the 
disordered and ordered phases. Above Tc, M^ = and MjJ = M 7^ 0. It is 
convenient to perform a rotation in the isospace: 
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In the rotated basis (T^^) = M and the Hamiltonian has the form: 



1 






m' 



(26) 



We then use the equations of motion for the operators T^ = T^ ± iT^: 
± if^ = 2t^T^ + i E y^^T': (t'^ + T-) , (27) 



in which we substitute T^ by its average value, M. This approximation, 
equivalent to the decoupling of the Green functions containing more than 
two isospin operators |3^, ^, gives a closed system of equations. 



uj-2t 



\/(q) ((T'+|r-))-fv^(q)((r-|r-)) 



2M 



(28) 



f v^(q)((r'+|r-)) + [uj + 2t± + f \/(q)] {{r-\r-)) = o 

for the retarded Green functions {{A\B)) = {{A\B))^^q, defined by 

POO 

((A|5)).,q = -I / t/te^(-+*^)*Ee"*''"'" ([A 5]) • 

JO n 

In this way we get: 
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(30) 
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Furthermore, 

(31) 

The Green functions have poles at (plus/minus) the energy of the isospin 
excitation with wave vector q: 



^q = V2tx {2t± + MV{q)). (32) 

The band of isospin excitations has its bottom at the wave vector q = Q, 
where V{q) has its minimum. 

Combining Eqs. (|30|) and ( pT] ) we obtain: 

((TIT^)).,, = ((T-IT-)) = ^\^^ . ^. (33) 

iu^ — £/q + tosignujj 

In the ordered phase below T^ we introduce on each site n an angle (^n by 

f M^ = MnCOS0n , . 

\ M^ = M„sin0„ ' ^"^ > 

where we used the original basis. Then the first of the self-consistency equa- 
tions Eq.dl) reads: 

2t_L sin 0n + cos 0n ^ Kim^^m siu 0^ = 0. (35) 

m 

We now chose the basis in the isospin space, such that the ^-axis is ori- 
ented along (Tn): 

T'l = sin^nT^ - coscpnT^, 

T'l = Ty, (36) 

T'n = cos (pnT^ + sin 4>nT^. 

Then (T'^) = M„ and (T'l) = (T'l) = 0. 

Using Eq.(p5|) the linearized equations of motion for T'^ can be written 
in the form: 

± zT^± = -^T't + ^M„ cos <f)n E ^nm COS <j)^ (t'+ + T'^) . (37) 

COS0n 2 ^ ^ ^ 
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Both for the zigzag and the chain ordering, 

0„ = e^Q-^>, (38) 

and \(f)n\ = <t>- Then, Eq.(p5D gives M^ = M and 

cos0|\/ (Q)| 

(c/. Eg ■ (|lO|) ) ■ Furthermore, from Eg . (1371) we see that the Green functions in 
the ordered phase can be obtained from those above Tc by substituting 2t±_ 
by 7^2 and 1^(q) by ^^(q) cos^ 0, e.g., 



fu; + ^ + f cos2 0V^(q) 

((T'+ir-)).,, = 2mA — r%2l-,- ^ (40) 



where the energy of the isospin excitations is given by 

2t, 



E. 



^ cos I 



.1--^^^. (41) 



Since below T^ < cos0 < 1, the isospin gap in the ordered phase, A = 
Eq = 2t±ta.n(f) > 0. The dispersion relations Eg . (13^) and Eg. (|4lD can be 
combined to give Eq.([T5|) valid both above and below Tc. 

Using Eg . (|35|) , one can check that the Green functions of T'^ and any 
other operator vanish. Then, ((T^|T^))(^^q = cos^ 0((T'^'|T'^')), so that below 



-CI 
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